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Abstract 


For ages flywheels have been used to achieve smooth operation of machines. The early models 
where purely mechanical consisting of only a stone wheel attached to an axle. Nowadays flywheels 
are complex constructions where energy is stored mechanically and transferred to and from the 
flywheel by an integrated motor/generator. The stone wheel has been replaced by a steel or composite 
rotor and magnetic bearings have been introduced. Today flywheels are used as supplementary UPS 
storage at several industries world over. Future applications span a wide range including electric 
vehicles, intermediate storage for renewable energy generation and direct grid applications from 
power quality issues to offering an alternative to strengthening transmission. 

One of the key issues for viable flywheel construction is a high overall efficiency, hence a reduction 
of the total losses. By increasing the voltage, current losses are decreased and otherwise necessary 
transformer steps become redundant. So far flywheels over 10 kV have not been constructed, mainly 
due to isolation problems associated with high voltage, but also because of limitations in the power 
electronics. Recent progress in semi-conductor technology enables faster switching and lower costs. 
The predominant part of prior studies have been directed towards optimising mechanical issues 
whereas the electro technical part now seem to show great potential for improvement. An overview of 
flywheel technology and previous projects are presented and moreover a 200 kW flywheel using high 
voltage technology is simulated. 
© 2006 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Several hundred years ago pure mechanical flywheels where used solely to keep machines 
running smoothly from cycle to cycle, thereby render possible the industrial revolution. During 
that time several shapes and designs where implemented, but it took until the early 20th 
century before flywheel rotor shapes and rotational stress were thoroughly analysed [1]. Later 
in the 1970s flywheel energy storage was proposed as a primary objective for electric vehicles 
and stationary power backup. At the same time fibre composite rotors where built, and in the 
1980s magnetic bearings started to appear [2]. Thus the potential for using flywheels as electric 
energy storage has long been established by extensive research. 

More recent improvements in material, magnetic bearings and power electronics make 
flywheels a competitive choice for a number of energy storage applications. The progress in 
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power electronics, IGBTs and FETs, makes it possible to operate flywheel at high power, 
with a power electronics unit comparable in size to the flywheel itself or smaller. The use of 
composite materials enables high rotational velocity with power density greater than that 
of chemical batteries. Magnetic bearings offer very low friction enabling low internal losses 
during long-term storage. High speed is desirable since the energy stored is proportional to 
the square of the speed but only linearly proportional to the mass. 

There are a number of attributes that make flywheels useful for applications where other 
storing units are now used. 


e High power density. 

e High energy density. 

e No capacity degradation, the lifetime of the flywheel is almost independent of the depth 
of the discharge and discharge cycle. It can operate equally well on shallow and on deep 
discharges. Optimizing e.g. battery design for load variations is difficult. 

e The state of charge can easily be measured, since it is given by the rotational velocity. 

e No periodic maintenance is required. 

e Short recharge time. 

è Scalable technology and universal localization. 

e Environmental friendly materials, low environmental impact. 


One of the major advantages of flywheels is the ability to handle high power levels. This 
is a desirable quality in e.g. a vehicle, where a large peak power is necessary during 
acceleration and, if electrical breaks are used, a large amount of power is generated for a 
short while when breaking, which implies a more efficient use of energy, resulting in lower 
fuel consumption. 

Individual flywheels are capable of storing up to 500 MJ and peak power ranges from 
kilowatts to gigawatts, with the higher powers aimed at pulsed power applications. 

The fast responstime in flywheels makes them suitable to balance the grid frequency. As 
the energy contribution from more irregular renewable energy sources increases, this can 
be an important quality which will grow in importance [3]. 

The tools used for motor/generator design are also continuously improving to yield a 
more correct picture of the induction process. Powerful computer programs, where full 
electromagnetic field calculations are considered, have reduced a number of limitations 
and approximations at the design stage. Along with technical progress, especially 
regarding high voltage generators, each machine can be designed to match the physical 
conditions of the energy source and of the load. In this way the electric efficiency of new 
machines can be increased significantly. 

High current yields substantial resistive power loss in the stator cables. To maximize the 
conductor area and thereby reduce the cable resistance, conventional generators use 
rectangular conductors. In theory it would be advantageous to build a generator that 
produces high voltage and low current, as the resistive power loss in the stator cables is 
proportional to the square of the current. Such a generator needs insulated circular 
conductors, for example conventional high voltage extruded solid dielectric cables [4]. This 
new class of generators is called Powerformer!™ [5-8]. 

The aim of this article is to give an overview of flywheel technology, its applications and 
present development. Furthermore the possibility of using high voltage motor/generators 
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in flywheels are discussed and a generator constructed with XLPE- suited for flywheel 
mounting is simulated and data are presented. 


2. Flywheel basics 
2.1. Energy storage in flywheels 


A flywheel stores energy in a rotating mass. Depending on the inertia and speed of the 
rotating mass, a given amount of kinetic energy is stored as rotational energy. The flywheel 
is placed inside a vacuum containment to eliminate friction-loss from the air and 
suspended by bearings for a stabile operation. Kinetic energy is transferred in and out of 
the flywheel with an electrical machine that can function either as a motor or generator 
depending on the load angle (phase angle). When acting as motor, electric energy supplied 
to the stator winding is converted to torque and applied to the rotor, causing it to spin 
faster and gain kinetic energy. In generator mode kinetic energy stored in the rotor applies 
a torque, which is converted to electric energy. Fig. 1 shows the basic layout of a flywheel 
energy storage system [9]. Apart from the flywheel additional power electronics is required 
to control the power in- and output, speed, frequency etc. 

The kinetic energy stored in a flywheel is proportional to the mass and to the square of 
its rotational speed according to Eq. (1). 

l 5 
Ex, = =I (1) 
2 
where Ex is kinetic energy stored in the flywheel, Zis moment of inertia and œ is the angular 
velocity of the flywheel. The moment of inertia for any object is a function of its shape and 
mass. For steel rotors the dominant shape is a solid cylinder giving the following 
expression for I: 


1 1 
I= z rm= 5 nap (2) 


Bearing 


Containment 


Flywheel rotor 


Motor/generator. Rotor 
Motor/generator. Stator 


Vaccum or Verry low 
preassure 


Bearing 


Fig. 1. Basic layout of a flywheel energy storage system [9]. 
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where r is the radius and a is the length of the cylinder, m represents the mass of the 
cylinder and p is the density of the cylinder material. The other dominating shape is a 
hollow circular cylinder, approximating a composite or steel rim attached to a shaft with a 
web, which leads to Eq. (3). 


1 1 
I= qin + r?) =] nap(r4 — r?) (3) 


Eq. (1) shows that the most efficient way to increase the stored energy is to speed up the 
flywheel. The speed limit is set by the stress developed within the wheel due to inertial 
loads, called tensile strength o. Lighter materials develop lower inertial loads at a given 
speed therefore composite materials, with low density and high tensile strength, is excellent 
for storing kinetic energy [10]. The maximum energy density with respect to volume and 
mass, respectively, is: 


ey =K; em =Ko/p (4) 


where e, and em is kinetic energy per unit volume or mass, respectively, K is the 
shapefactor, o is maximum stress in the flywheel and p is mass density. In case of planar 
stress, if the height of the disk is small compared with the diameter, and a homogenous 
isotropic material with Poisson ratio of 0.3, i.e. steel, is used, the K factors are given in 
Table 1 [11]. 

In a three-dimensional object there will be three-dimensional interaction of material 
stresses. For a rotor constructed with a non-isotropic material, like fibre-reinforced 
composite, that stress interaction will limit the practical dimensions possible. Taking into 
account safety issues the resultant flywheel design is based on a hollow cylinder, in which 
material stresses created by three-dimensional effects are minimized. In short designs, the 


Table 1 
Shape-factor K for different planar stress geometries 


Fly wheel geometry Cross section Shape factor K 


Disc Lif 1.000 


Modified constant stress disc AY 0.931 


Conical disc ZOZ 0.806 
Flat unpierced disc ezzzzzzz TTI TA 0.606 


Thin firm as aT 0.500 


Shaped bar ee i 0.500 
Rim with web Z; 0.400 


Single bar | 0.333 


Flat pierced bar ZZZ 7777777777 0.305 
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two stresses of primary concern are the radial stress and the hoop stress, Fig. 2. For an 
isotropic material the radial stress is expressed by Eq. (5). 


3+0 rere 
of 5 2 
Fr =~ pw ah = at 


(5) 


where p is the mass density, œ is the rotor speed, v represent the Poisson ratio rọ is the outer 


radius of the rotor, r; is the inner radius of the rotor and r represent any radius within the 
rotor. 


The hoop stress is expressed by Eq. (6) 


n= (Re - 4?) 


5 © 


Table 2 presents characteristics for common rotor materials [12]. 

By making fiber reinforced composite rotors with circumferentially oriented fibers, the 
flywheel is more likely to develop circumferential cracks, which are much less likely to 
produce free-flying projectile fragments in case of a catastrophic failure. In most designs a 
rotational speed drop of 50% is allowed, thus the available energy is 75% of the stored 
energy, in other words the depth of discharge is 75%. Overall the flywheel geometry and 


speed determines the energy storage capability, whilst the motor/generator and power 
electronics determines the power capabilities. 


\ \ gadd 
-3 Æ Al ess 


Ws 
3s 


Fig. 2. Radial- and hoop stress in a short hollow cylinder rotating about its axis with angular velocity œ. 


Table 2 

Data for different rotor materials 

Material Density Tensile strength Max energy density (for | kg) Cost ($/kg) 
(kg/m*) (MPa) 

Monolithic material 7700 1520 0.19 MJ/kg = 0.05 kWh/kg 1 

4340 Steel 

Composites 

E-glass 2000 100 0.05 MJ/kg = 0.014 kWh/kg 11.0 

S2-glass 1920 1470 0.76 MJ/kg = 0.21 kWh/kg 24.6 

Carbon T1000 1520 1950 1.28 MJ/kg = 0.35 kWh/kg 101.8 

Carbon AS4C 1510 1650 


1.1 MJ/kg = 0.30kWh/kg 31.3 
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2.2. Magnetic bearings 


Mechanical bearings used in the past cannot, due to the high friction and short life, be 
adapted to modern high-speed flywheels. Instead a permanent or electro permanent 
magnetic bearing system is utilized. Electro permanent magnetic bearings do not have any 
contact with the shaft, has no moving parts, experience little wear and require no 
lubrication. It consists of permanent magnets, which support the weight of the flywheel by 
repelling forces, and electromagnets are used to stabilize the flywheel, although it requires 
a complex guiding system. An easier way to stabilize is to use mechanical bearings at the 
end of the flywheel axle, possible since the permanent magnet levitates the flywheel and, 
thus, reduce the friction [13,14]. The best performing bearing is the high-temperature 
super-conducting (HTS) magnetic bearing, which can situate the flywheel automatically 
without need of electricity or positioning control system. However, HTS magnets require 
cryogenic cooling by liquid nitrogen [12]. 


3. Flywheel technical considerations 


For decades, most engineers have used the concept of storing kinetic energy in a 
spinning mass to smooth their operation. Until recently the vast majority constituted of 
steel wheels coupled with a motor/generator, where the high rotary inertia allowed long 
ride-through time without significant decrease in flywheel rotational speed. Since the 
change in rotational speed directly reflects the electrical frequency the power delivery of 
those flywheels rarely exceeded 5% of the stored energy. 


3.1. Motor/generator 


Requirements for standardized electric power have made most flywheel system designers 
elect variable speed AC generators (to accommodate the gradual slowing of the flywheel 
during discharge) and diodes to deliver DC electricity. The two major types of machines 
used are the axial-flux- and the radial-flux permanent magnet machines (AFPM and 
RFPM, respectively). There are numerous alternatives for the design of an AFPM machine 
such as internal rotor, internal stator, multidisc, slotted or slot-less stator, rotors with 
interior or surface-mounted magnets [15-18]. Unlike radial machines, axial machines can 
have two working surfaces. Either two rotors combined with one stator or one rotor 
combined with two stators. The benefit of using a two surface working machine is the 
increase in power output [17,19]. The axial machines seem to have more advantages over 
the radial such as, a planar adjustable air gap and easy cooling arrangements, which is 
important when working under low-pressure conditions [20]. Fig. 3a shows a one-rotor 
two stator AFPM configuration without the cable winding in the stators. It can be seen 
that the permanent magnets are an integral part of the flywheel rotor and the stators are 
fixed to the housing. 

A permanent-magnet axial-flux motor, resembling the 200kW motor/generator 
simulated later in this article, is analysed by Furlani [21]. An analytical expression for a 
two dimensional field solution of the magnetic field is presented for a simple pole structure. 
The permeability for the magnetic field return path is set to infinity. 

Much attention has been directed towards optimising radial gap machines [22-24]. In a 
RFPM machine the magnets can be surface mounted on the rotor axle surrounded by the 
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Fig. 3. (a) show an AFPM machine arrangement and (b) show an RFPM machine arrangement. 


Fig. 4. Cross section of internal dipole array for n = 8. 


stator, as in Fig. 3b, or mounted in a ring enclosing the stator. The radial flux machine is 
mostly used in small-scale high-speed machines, where the tensile strength of the 
permanent magnets demands placing close to the rotating axle. 

Another type of motor/generator is the internal-dipole, Halbach-type magnet arrays. 
Where the PM array rotates with the flywheel and interacts with a set of stationary coils to 
produce torque. In a Halbach array, n PM segments forming a cylindrical shell about an 
axis create the internal dipole. A cross section of an internal dipole array with 8 segments, 
where M is the magnetization is shown in Fig. 4, inside a single turn two-phase stator is 
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also shown. The Halbach type motor can also be of multi-pole type. One of the advantages 
of this configuration is the low external magnetic field produced when a steel rim is placed 
outside the magnets. Within the shell a dipole configuration creates a uniform flux B [25]. 


B = Bem log 2 K (7) 
rı 
sin (2) 
=e (8) 
p 


where B is the resultant uniform flux, B,em is the remanent flux in permanent magnets, r2 
equals outer radius and rı equals inner radius and p is the number of poles. 

The torque generated is then proportional to the ampere turns and radius from the 
rotational center [25]. 


T = Bilr (9) 


where T equals torque, B is the magnetic flux, / is the length of the conductor, i is the 
current and r equals the radius. 
Machines with operating voltage in the range of 70-400 V have been built [26,27]. 


3.2. High voltage 


Even though different kinds of flywheels constructed today benefit from the recent 
progress in technology, there is one thing all of them have in common, the inability to 
directly produce high voltage (>36kV). So-called ‘high voltage’ flywheels have been 
constructed. However, the highest voltage attained so far is a 10-pole permanent magnet 
machine with a continuous voltage of 6.7kV and a peak voltage of 10kV constructed in 
2001 [28]. The result is that for true high voltage applications a transformer has to be used, 
introducing more unwanted losses. 

Apart from the PM motor/generator used in almost all flywheels there is also the 
possibility of using a Synchronous Reluctance Motor/Generator. In 1996 a 60kW 
flywheel, utilizing this motor, was developed [29]. Table 3 shows advantages/disadvantages 
with PM and induction machines. 


3.3. Number of poles 


The choice of number of poles to be used in a machine is essential to the overall 
performance. Two pole motor/generators are most common in high-speed machines, 
mainly to keep the voltage down but it also has other good properties. Depending on 
axialor radial flux configuration a multi-pole rotor can experience substantial electro- 
magnetic axial or radial forces generated by the stator winding, if there is a net attractive 
force between a pole-pair and the stator. In a two-pole rotor, however, the only two poles 
are directly opposite one another resulting in a net force on the rotor of approximately 
zero. Eliminating these forces reduces the load requirements on the bearings, which is 
particularly important if magnetic bearings are used. 
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Table 3 


Advantages and disadvantages with PM and induction machines 


High voltage PM machine 


PM machine 


B. Bolund et al. / Renewable and Sustainable Energy Reviews 11 (2007) 235-258 


Induction machine 


Advantages + High overload capability, + Magnetic field is + No concern with 
due to low load angle and produced without excitation demagnetization 
low stator current losses 
+ Magnetic field is + Less complex rotor + No excitation field at zero 
produced without excitation design, no need of electric torque, hence no 
losses wires in the rotor electromagnetic spinning 
losses 
+ Less complex rotor + Possible to achieve a + Can be constructed from 
design, no need of electric higher overall efficiency high-strength low-cost 
wires in the rotor materials 
+ Possible to achieve a 
higher overall efficiency 
Disadvantages — Risk of demagnetization — Risk of demagnetization — RP, transformer and 


and a decreasing intrinsic 
coercivity with increasing 
temperature 

— Machines with iron in the 
stator experience 
electromagnetic spinning 
losses at zero-torque 

— The low tensile strength of 
PM materials require 
structural support against 
centrifugal forces, leaving 
constraints on the design of 
high-speed, high-power 
rotors 


and a decreasing intrinsic 
coercivity with increasing 
temperature 

— Machines with iron in the 
stator experience 
electromagnetic spinning 
losses at zero-torque 

— The low tensile strength of 
PM materials require 
structural support against 
centrifugal forces, leaving 
constraints on the design of 
high-speed, high-power 
rotors 


rectifying losses in the 
electromagnets during field 
excitation 

— More complex rotor 
design, due to the need of 
wires and electric 
brushconnection to the rotor 
— Rotor brushes require 
maintenance 


— Poor overload capability 
due to the high stator current 


3.4. Power electronics 


A brushless permanent magnet generator (in a flywheel) produces variable frequency AC 
current. In most applications though, the load requires a constant frequency making it 
necessary to first rectify the current and then convert it back to AC. Power converters for 
energy storage systems are based on SCR, GTO or IGBT switches. In an early stage of 
energy storage utility development, SCRs where the most mature and least expensive 
semiconductor suitable for power conversion. SCRs can handle voltages up to 5kV, 
currents up to 3000A and switching frequencies up to 500Hz. Due to the need of an 
energized power line to provide the external on/off signal to those switches they where 
replaced with GTOs, which do not depend on an energized line to function. The GTO 
device can handle voltages up to 6kV, currents up to 2000 A and switching frequencies up 
to 1 kHz. In the last several years IGBTs has emerged, Fig. 5. The IGBT is a solid-state 
switch device with ability to handle voltages up to 6.7kV, currents up to 1.2kA and most 
important high switching frequencies [30,31]. 
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Collector 


— 


=... Emitter sess 
Gate Gate 


Fig. 5. IGBT schematic. 


The technique used to produce AC current from DC is called Pulse-Width Modulation 
(PWM). Pulses of different length are applied to the IGBTs in the inverter, causing the DC 
current to be delayed by the inductive load and a sine wave is modulated [32]. A fast 
switching frequency in the power converter improves emulation of a sine wave mainly by 
eliminating some of the higher order harmonics. To reduce the harmonic content even 
further a filter, consisting of capacitors and inductors can be connected on the AC side of 
the output. 

Electromechanical energy conversion systems have been explored in order to make them 
more fault-tolerant [33]. By using three single-phase inverters instead of one compact setup 
a more flexible design is achieved along with advantages like the ability to operate even in 
the event of a single-phase fault. 


3.5. Work done to date 


3.5.1. Small-scale 

Small-scale flywheel energy storage systems have relatively low specific energy figures 
once volume and weight of containment is comprised. But the high specific power possible, 
constrained only by the electrical machine and the power converter interface, makes this 
technology more suited for buffer storage applications. Development of alternative dual 
power source electric vehicle systems that combine a flywheel peak power buffer with a 
battery energy source has been undertaken [18,34]. 


3.5.2. Peak power buffers 

The uses of a flywheel as power buffer in an electric vehicle can significantly reduce the 
peak currents drawn from the ordinary storing supply e.g. battery. Elimination of the peak 
currents will prolong the battery life [34]. 


3.5.3. Wind-diesel generator with a flywheel energy storage system 

In the year 2000 a simulation of a Wind-Diesel generation plant together with a kinetic 
energy storage unit was presented and the construction of it was undertaken. The goal of 
this system is a unit where the regular wind oscillations are compensated by the diesel 
generator and the flywheel. The 0.6kWh, 50kW flywheel is able to supply active and 
reactive power to compensate both frequency and voltage of the network. The unit is 
designed to supply total power during a period of 1.8 min with a rated voltage 750 V and a 
maximum current of 102 A [35]. 
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3.5.4. Flywheel for photovoltaic system 

A doubly salient permanent magnet (DSPM) motor flywheel energy storage for building 
integrated photovoltaic (BIPV) system was simulated in 2001. By adding a flywheel to a 
BIPV equipped building situated in Hong Kong, the load supply time can be prolonged 
from 9 a.m. to 3 p.m. to 8 a.m.-beyond 6 p.m. [36]. 


3.5.5. Harmonics 

Different flywheel systems for compensating harmonics in low voltage (~400 V) power 
networks have been compared and analysed. Up to the eleventh harmonic a decrease of 
about 50% was accomplished [37]. 


3.5.6. Flywheel in distribution network 

A 10MJ flywheel energy storage system, used to maintain high quality electric power 
and guarantee a reliable power supply from the distribution network, was tested in the year 
2000. The FES was able to keep the voltage in the distribution network within 98—102% 
and had the capability of supplying 10kW of power for 15 min [38]. 


3.5.7. High power UPS system 

A 50MW/650 MJ storage, based on 25 industry established flywheels, was investigated 
in 2001. Possible applications are energy supply for plasma experiments, accelerations of 
heavy masses (aircraft catapults on aircraft carriers, pre-acceleration of spacecraft) and 
large UPS systems. The 50 MW peak power can be supplied for about 13 s, with an overall 
efficiency of 91-95%. The flywheels are connected in parallel to a 1200 V DC-link. Similar 
PM flywheels have previously been tested in urban traffic busses and rail systems with a 
resulting energy save of up to 40% [39]. 


3.5.8. UPS system 

A Case study on an existing medium voltage network has been carried out, in which 
different disturbance scenarios have been simulated (voltage dips, start-up etc). The idea 
was to connect four 1.6 MVA flywheel based dynamic UPS systems combined with a diesel 
generator to the 20kV distribution network, thereby improving the power quality. The 
simulation results indicate that the approach is feasible, and show a significant 
improvement in power quality. Typically, the voltage dips is divided by a factor 3 even 
in the worst cases. A transformer is required between the flywheel storage system and the 
medium voltage network [40]. 


3.5.9. Aerospace applications 

A two pole, three-phase PM synchronous motor/generator coupled to a flywheel have 
been simulated. The flywheel storage unit is intended to replace a battery storage unit 
onboard the International Space Station. The motor is rated to 7kVA, 80 V and 50A and 
1000 Hz. A comparison between flywheel and NiH, battery systems for an EOS-AMI type 
spacecraft has shown that a flywheel system would be 35% lighter and 55% smaller in 
volume [41]. 


3.5.10. High voltage stator 
A 10 pole PM machine with a continuous voltage of 6.7kV and peak voltage of 10kV 
was constructed in 2001, for use in hybrid electric combat systems. The 25 MJ flywheel is to 
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produce a continuous power of 350 kW (loss 2.4kW) as well as intermittent 5 MW pulses, 
the idling loss is calculated to be around 250 W. To handle the heat produced, the stator is 
cooled by 70 or 90°C oil. The insulation on the cables constitute of filler epoxy (type 1) and 
FEP tubes surrounded by epoxy (type 2) see Fig. 6. The overall size and weight of the 
system is 0.28 m° and 519 kg [28]. 


3.5.11. Of the shelf systems 
The operating voltages of three available flywheel systems from different companies can 
be found in Table 4. 


3.6. Losses 


The integrated motor/generator is usually of a rotating-field deign, where the field is 
supplied either by electromagnets or by rare-earth permanent magnets. The properties of 
high field permanent magnets yield flux densities high enough to enable machines with 
slotless armature windings, also known as air-gap windings, without a magnetic stator core 
[42]. Absence of a ferromagnetic material in the stator has two major impacts on the 
performance of a motor/generator. First of all the low permeability will quickly reduce the 
magnetic field strength when moving away from the magnet. As a result the induced 
voltage becomes lower, thus also diminishing the generated power. Second, there will not 
be any heat loss in the stator core due to hysteresis effects. The hysteresis loss is otherwise 


Feb oil flow tubes 


Filler Epoxy 


Cable with internal Epoxy 
0.020 FEP Roll Cover 


Fig. 6. Cross section of the 6.7kV Stator [28]. 


Table 4 

Operating voltage of three different flywheel systems [12] 

Developed by: Type Power (kW) Voltage (V) 
EMAFER Medium power 300 < 1000 
Magnet motor Medium power 150 300-800 
Magnet motor High power 5000 — 


URENCO Medium power 100 600-750 
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always present when exposing a ferromagnetic material to magnetic flux [43]. It is clear 
that higher frequencies convey more losses and that hysteresis loss in the stator core will 
have severe impact during long time (stand-by) energy storage in a flywheel. Without 
hysteresis loss the stand-by losses are very small and limited to those of leak eddy currents 
and bearing losses. 

The low tensile strength of the magnets compared to that of the composite 
flywheel limits their placing to the vicinity of the hub. As a consequence the number of 
poles, and therefore the rate of change of magnetic flux, must be carefully selected in 
order to achieve the desired voltage. Table 5 shows the tensile strength for common 
magnetic materials [44]. 

Traditionally used ferrites do not, due to their low conductivity, give rise to induced 
eddy currents on the surface. Some of the sintered rare earth materials, however, have large 
conductivity and therefore suffer from such problems. Eddy currents on the surface of the 
magnets arise when the magnetic field from the stator interacts with the magnets and are 
roughly: 


J _ T AP sator 
magnet 5S7 jy 
Imagnet dt 


(10) 
where Jmagnet is the surface current density on the magnet, o is the conductivity of the 
magnet, /magnet iS the length of the magnet, and Pstator is the magnetic flux from the stator. 
Since the magnetic flux form the stator is proportional to the current going through the 
stator windings and eddy current losses are squarely dependent on frequency, it is 
necessary to minimize the current harmonics. 

A great deal of the total losses from the motor is ohmic (I R) losses in the stator 
winding. It is clear that those can be diminished either by increasing the amount of 
conducting material (usually copper), thereby decreasing the resistance, or by decreasing 
the current in the stator. There are obvious drawbacks associated with increasing the 
amount of conducting material such as increasing weight, cost and space. Induced high 
frequency eddy currents in the stator may also increase depending on the configuration. 
Decreasing the current in the stator inevitably leads to higher voltage if the overall power is 
to be maintained. So far it has been impossible to increase the voltage due to the risk of an 
electric breakdown, but if a higher voltage can be handled the copper losses can be 
decreased. 


Table 5 
Data for different magnetic materials [44] 


Material Density Tensile strength Remanence 
(kg/m?) (MPa) (T) 
Sintered Neodymium—Iron—Boron (Nd—Fe-B) 7400-7600 80 1.08-1.36 
Sintered Samarium-Cobalt 8000-8500 60 0.75-1.2 
Sintered Ferrite 4800-5000 9 0.2—0.43 
Injection molded composite (Nd-Fe-B) 4200-5630 35-59 0.40-0.67 
Compression molded composite (Nd-Fe-B) 6000 40 0.63-0.69 


Injection molded composite Ferrite 2420-3840 39-78 0.07-0.30 
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Fig. 7. Simple illustration of how to mount a flywheel in a vehicle for minimizing the effect of gyroscopic torques. 


3.7. External gyroscopic aspects 


For flywheels situated in a vehicle, satellite or space station the gyroscopic forces are 
important. In a satellite it is possible to make use of the gyroscopic force from a flywheel 
by letting it provide torque to the spacecraft for altitude control, thereby diminish the 
overall weight and increase the effective use of energy. In space stations and vehicles the 
flywheel batteries are controlled as pairs and situated to rotate in opposite direction to not 
produce any net torque [45,46]. Another way to cope with the interaction of gyroscopic 
forces in a vehicle when using just one flywheel is to place the flywheel in a gimbal system, 
thereby eliminating most of the gyroscopic torque. The gimbal system works in the same 
way as a cup-holder does in a vehicle, which means that the vehicle can turn and lean 
without tilting and twisting the position of the flywheel [47]. Fig. 7 shows a simplified 
picture of the gimbal system. 


3.8. Safety 


Some sort of inertial containment system becomes necessary to minimize the collateral 
damage from a failed flywheel. Reasons for failure could be, crack growth from material 
flaws created during manufacture, bearing failure or external shock loads. For large 
flywheels the vacuum chamber acts as a first safety enclosure in a multiple-barrier 
containment system to prevent rotor debris from flying free. The next barrier system design 
can include thick steel, concrete chambers and/or underground vaults. Small portable 
flywheels cannot utilize bulky containments like an underground vault. Instead the rotor is 
designed to fail safely (described above) in which case a vacuum vessel provides sufficient 
protection. It is also possible to simply place the FES units in restricted areas similar to 
what is done with conventional turbines that operates in electric power plants. Most 
machines have a vertical rotation axle, but horizontal machines also occur [48]. A vertical 
axle minimizes the possibility of mass centre displacement witch can lead to instabilities 
and damage the flywheel. 


4. Simulation of motor/generator 
4.1. Mathematical model of the generator 


When a generator is simulated a two-dimensional model of the generator cross-section 
geometry is created. The geometry is based on straight lines and circular arcs. The 
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geometric domains are assigned a material with corresponding material properties such as 
resistivity, permeability, coercivity, sheet thickness, price etc. Voltage, current and thermal 
sources are given as scalars or by circuit equations. 

A two-dimensional finite element method (FEM) is used in the calculations. To ensure 
high accuracy and fast computations the mesh is made more detailed in the air-gap and in 
the stator teeth and coarser in the yoke and the rotor rim. The accuracy can be set to 
different levels. To account for three-dimensional effects, coil-end reactances and 
resistances are calculated. 

When calculating the induction in the stator the displacement current, 0D/0t, can, due 
to its long wavelength, be neglected. The displacement current is also directed in radial 
direction in the insulating dielectric material surrounding the stator cables and will thereby 
not contribute to any induction. 

Without the displacement current Ampere’s law can be written as: 


Vx H=j (11) 


where H is the magnetizing field and j represent the free current density. Material 
behaviour is described by Eqs. (12) and (13). 


B = Uy MoH (12) 
where uuo is the magnetic permeability and B is the magnetic flux, Ohm’s law, 
J=oE (13) 


where ø is the conductivity and E is the electric field. The B-field can be expressed by the 
vector potential, A, as: 


B=VxA (14) 
Combining Eqs. (11)-(14) gives: 


vx (= Vx A) =cE (15) 
LU; Ho 
Faradays law: 
oB 
Vx E = —-— 1 
x F (16) 
in combination with Eq. (14) and Helmholtz’ theorem gives: 
0A 
E=-—-VWV 1 
a (17) 


where V is the scalar potential. The two-dimensional nature of the model enables the 
vector potential to be expressed as: 


A= A,r, 9, À (18) 
which together with Eqs. (15) and (17) results in: 
0A, 1 OV 


The term 0V/0z is traditionally called applied potential and can be given as an 
inparameter corresponding to a current density in z-direction. The time derivative of the 
vector potential in Eq. (19) is related to the penetration of a magnetic field in a material, 
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also called the skin effect. The skin depth in a material, 6,4in, is given by: 
1 


Ôskin X —— (20) 
where f is the electric frequency of the generator. 
A current source representation of Eq. (15) can be written as 
1 
Vx-VxA=J+Jm (21) 
u 


where J is the source current density given by Eq. (11) and the current density Jm is given 
by Eq. (23) [49]. In the simulation tool Eq. (21) is solved for a number of different discrete 
rotor positions under transient conditions, where the distance between those positions is 
given by the rotor speed. The rotor and stator are connected with varying boundary 
conditions. Symmetries in both geometry and electromagnetic field enables the generator 
to be presented by a two dimensional unit-cell with periodical boundary conditions. 
Depending on the stator slot pitch the unit cell can include one or more rotor poles. In the 
case of a permanently magnetized rotor, as in this case the length of the permanent magnet 
is determined by iteration to give a sufficient induction. 
Permanent magnets are modelled by a surface current density, Jm,s, determined by: 


Jms =AxM (22) 


where M is the magnetization. 

If stator steel is used there will be distributions of non-linear magnetization and current 
in the stator. In that case the magnetization current, Jm, in the volume representing the 
stator steel is given by: 


JIm=VxM (23) 


The Flywheel motor/generator simulated in this article uses no stator steel and therefore 
no attention needs to be paid to non-linear currents and magnetization in the stator. 
Thermal distribution in the generator is determined by Fick’s law and the continuity 
equation for heat: 


oT 
kWT = Ot = —S heat (24) 


where k is the diffusion coefficient and Jheat is a heat source i.e. ohmic current losses, 
magnetic losses and external cooling. In addition to the field equations appropriate initial, 
boundary and jump conditions should be added. The values in the magnetization curves, 
BH-curves, for all materials have been experimentally derived by the Epstein method. 

Electromagnetic losses in the generator consist mainly of ohmic and eddy current losses 
in the stator copper winding. The mesh in the FEM-solver has higher mesh density in areas 
of special interest such as in the air gap and in the stator. Among the calculation results are 
magnetic field plots, temperature distributions, iron and copper losses, load angles and 
reactances. 


4.2. Assumptions and design objectives 


A 200 kW three-phase PM generator intended for use in a flywheel storage unit, situated 
in for example a bus, has been simulated. To minimize stand-by losses an ironless stator 
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was preferred. To achieve high power density in a generator with an ironless stator, an 
axial flux topology where the stator cables are situated between one magnet and one 
ferromagnetic steel rim was chosen. The magnets used are of neodymium—iron-boron, 
NdFeB, type. As this material is conductive it is necessary to limit the current 
harmonics from the stator. Otherwise induced surface eddy currents in the magnets 
will cause excessive loss and heat development, which could deteriorate magnet 
characteristics. Conventional high voltage extruded solid dielectric cables are used for 
winding the stator [4]. 

Models for the electric circuit are added to the electromagnetic and thermal field 
equations. The equations describing the stator (armature) circuit are, 


ia + ip +i, = 0 


Ua + Rsig + ise a up — Rsip i a = Va 
Ole Oi 
Ue + Ric F Le ~ Up Rsip ge = Veb (25) 


where ua, Up are the phase voltages, i,, i, are the phase currents, p is the coil end 
reactance (obtained from an separate three dimensional computation), and Vab, Vep are the 
line voltages of the external load circuit. The external load may be determined in more 
detail, by for instance including equations for rectifiers. However, a simple circuit with a 
resistive load is simulated here. 

The internal (phase currents, phase voltages, rotor and stator equations) and the 
external (outer load) circuit equations are added to the system. Sources are des- 
cribed by circuit equations. A symmetric three-phase winding of the stator is modeled. 
The phases are mutually separated by 27/3 electrical radians. The main advantage 
of a symmetric threephase system is that the output power is only slowly varying 
with time (i.e. the variation of rotation speed) provided a proper design is made for the 
generator. 


4.3. Motor/generator design 


An air-wound axial flux machine is modeled. The upper rim rotates with the same speed 
as the magnet and the lower rim. Fig. 8 shows a cross-section of one pole. 

Table 6 gives the properties of the generator. 

The machine is studied for a rotational speed of 8000rpm which corresponds to a 
tipspeed of 420 m/s for a carbon composite rotor with a diameter of 1 m. The machine has 
12 poles and is designed to fit inside the carbon composite cylindrical rotor. 

For simplicity the power factor cos@ is chosen to 1 in the simulations (in a traditional 
rotating machine it is typically around 0.8). In a real construction, the power factor could 
be optimized to the predominant rotational frequency. However, if the machine is designed 
to work at low load angle, as in this case, the power factor will be close to unity for the 
entire operating range. 

No assumption on the size or genometry of the carbon composite rotor, other 
than rotational frequency, are made. The machine is assumed to operate with a rotor 
that can store in the range of 5kWh which in practice would correspond to a rotor 
mass of 30 kg. 
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Steel 


Fig. 8. Cross-section view of one generator pole. 


Table 6 

Generator properties 

Power 200 kW 
Voltage 1kV 
Current 115.5A 
Power factor, cos b 1 

Load angle 4.3° 
Outer diameter 0.689 m 
Rpm 8000 
Max electric freq 800 Hz 


0 = 004 008 2 16 


am 


Fig. 9. Magnetic field in the air-gap for one pole. The grey line shows the B-field along a line in the middle of the 
air-gap and the black line is the fundamental of B. 


5. Results from simulation 


A plot of the magnetic field distribution along a line in the middle of the air-gap is shown 
in Fig. 9. The ripple causes a rapidly varying electromagnetic field on the rotor. This ripple 
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can cause vibration but also induce eddy currents which will contribute to losses as 
discussed above. Fluctuations can also be observed in the voltage output, Fig. 10. 

A total amount of 4.9 dm? of magnets are used. Fig. 11 shows the magnetic field in the 
generator under load condition. The losses are given in Table 7 and the voltage harmonics 
produced are given in Table 8. 
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Fig. 10. Stator voltage output from the three phases in generating mode at full load. 
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Fig. 11. Magnetic field distribution under full load in generator mode. 


B. Bolund et al. / Renewable and Sustainable Energy Reviews 11 (2007) 235-258 255 


Table 7 
Electric losses in the generator 


Losses (kW) 


Ohmic copper loss in core 0.1 
Ohmic end winding loss 0.29 
Eddy current loss in windings 0.79 
Table 8 


Voltage harmonics in the stator 


Phase harmonic I(%) 

1 100 

3 6.8 

5 0.4 
15 1.0 
17 1.9 
19 1.7 
21 0.6 
33 0.4 


The thermal computations are not yet complete, as some thermal insulation material 
properties outside the stator are not included. However, the maximum temperature of the 
generator is so far below 60° in all simulations. The increase in thermal load due to the low 
air-pressure around the motor/generator will be compensated by a water or air cooling 
system. A well-designed generator with low operational temperature opens for large 
overload capacities without operational worry of fast insulation degradation. The 
conductor area has been chosen to 16mm/”, which is sufficient to avoid strong Ohmic 
heating. 


6. Discussions 


The calculations are based on a two-dimensional electromagnetic and thermal field 
model. Three-dimensional effects are not fully taken into account, but the coil-end 
reactance is included. The generator physics is described by field equations, and a 
conventional simplified circuit description (‘equivalent circuit’) is thus not the basis for the 
simulations. 

One of the key issues for a working flywheel is to keep the induced eddy-currents to a 
minimum. Temperature rise in the rotor magnets due to eddy-currents from the phase 
harmonics needs to be analyzed. Different stator layouts will give rise to different amounts 
of harmonics, it is therefore crucial to simultaneously analyze eddy-currents and cooling. 
Although simulations show that the configuration used in the simulations work, there is 
still the need of simultaneous analyzing of the induced eddy-currents in the magnets as well 
as the need for experimental validation of the results. 

The voltage harmonics listed in Table 8 can be observed as the voltage ripple in Fig. 10. 
The stator current will contain about the same harmonics and since the induced eddy 
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currents in the rotor magnets will depend strongly on the shape of that stator current, it is 
necessary to try to minimize harmonics. 

For motor/generators operating under a low surrounding air-pressure, a high electric 
efficiency is crucial to avoid overheating. A high voltage ensures low currents which leads 
to low losses and opens for a large overload capability, up to 500 kW for this motor/ 
generator. The stator winding arrangement and rotor pole shape will have large influence 
on the overall performance. A very interesting choice is to have magnets mounted at both 
sides of the stator, which would strongly enhance the magnetic field in the air-gap. 


7. Conclusions 


Flywheel storage systems have been used for a long time. Material and semiconductor 
development are offering new possibilities and applications previously impossible for 
flywheels. The fast rotation of flywheel rotors is suitable for direct generation of high 
voltage. Thus for flywheel applications, the motor/generator part has a large upgrade 
potential. In this article a 200kW permanent magnet air gap winding motor/generator 
with axial flux has been simulated. This motor/generator setup incorporates high voltage 
technology with the use of NdFeB permanent magnets and an air wound stator. The 
simulation topology used was partly chosen for simulation simplicity and can most likely 
be enhanced. The simulated motor/generator is intended for a flywheel storage system 
situated in e.g. a bus. However, this flywheel technology is scalable and larger machines 
can be constructed for the applications of e.g. stabilizing the electric grid. 
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